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We investigate the possible effects of reionization by active sources on the cosmic microwave 
background. We concentrate on the sources themselves as the origin of reionization, rather than 
early object formation, introducing an extra period of heating motivated by the active character of 
the perturbations. Using reasonable parameters, this leads to four possibilities depending on the 
time and duration of the energy input: delayed last scattering, double last scattering, shifted last 
scattering and total reionization. We show that these possibilities are only very weakly constrained 
by the limits on spectral distortions from the COBE FIRAS measurements. We illustrate the 
effects of these reionization possibilities on the angular power spectrum of temperature anisotropies 
and polarization for simple passive isocurvature models and simple coherent sources, observing the 
difference between passive and active models. Finally, we comment on the implications of this work 
for more realistic active sources, such as causal white noise and topological defect models. We show 
for these models that non-standard ionization histories can shift the peak in the CMB power to 
larger angular scales. 


I. INTRODUCTION 

It is well known that the angular power spectra of temperature anisotropies and polarization of the cosmic microwave 
background (CMB) depend sensitively on the thermal history of the universe 01- Around the time of recombination 
(z ss 1100), when protons and electrons recombine into neutral hydrogen, the microscopic physical processes at work 
are relatively well understood and the calculation of the photon visibility function, which feeds into the angular power 
spectra, is relatively simple, at least in theories with passive fluctuations, such as inflation. 

However, even in these theories the universe must have become reionized, since there is no Lyman-a trough in 
distant quasar spectra (the Gunn-Peterson test ||). It is thought that this must be due to virialized objects, such 
as protogalaxies, massive stars and quasars, which formed relatively early in the history of the universe jl(j |l2j]. The 
microphysics of such processes is less well understood, but photoionization due to radiative objects can not happen 
earlier than when these objects have been created, which is believed to be at redshifts below z = 100, and it is 
actually thought more likely to have happened much later, after z ~ 30 |u|. If this is the case, then the actual 

observed CMB anisotropies and polarization will be a small perturbation on those calculated using the standard 
thermal history, just including recombination, due to the small optical depth of the time of reionization. Nonetheless, 
there are some potentially observable effects, particularly in the polarization |2(} |. 

As always the situation is much less clear in the case of actively generated perturbations, such as those from 
topological defects J2l] |23|]. In such models, one is forced to try to model highly non-linear processes from the time 
of defect formation to the present day, which is approximately 25 orders of the magnitude in expansion. Even with 
the most powerful super-computers available at present, it is difficult to achieve much more than a factor of 1000 
in expansion and hence extrapolations are necessary. Notwithstanding these difficulties there does seem to be a 
consensus at the present as to the predictions of the simplest models, cosmic strings and textures |24| [jl|], using 

the standard thermal history. It appears that flat universe models with critical matter density (fl m = 1) would require 
unacceptably large biases (~ 5) on 100h -1 Mpc scales to be consistent with the observed galaxy distribution [|d| p6| , 
although more exotic defect models may not have this problem |2^ [2^|. More acceptable models can be constructed 
in an open universe or one dominated by a cosmological constant ||2S|j3l| , |3^| . Even, if these models were to be ruled 
out by future observations, it is still important to investigate the possibility of more general active sources as the only 
credible alternative to inflation. 

The purpose of this paper is to investigate the possible effects of the active character of such sources for structure 
formation on the thermal history of the universe. The basic conceptual difference between active and passive models, 
such as inflation, is that the sources are present in the pre- and post-recombination plasma. This can create potentially 
large, local, non-linear perturbations, which can accelerate matter causing shock heating up to temperatures of a few 
million Kelvin in the baryons |3fj|| . In order to perform a convincing quantitative treatment of these effects one would 
have to incorporate the sources accurately into a full hydrodynamic simulation which has sufficient resolution to 
accurately model both the activation of ionization by the gravitational effect of the sources on the smallest scales 
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and also large enough to model the effects of an expanding universe. Obviously, the amount of computer resources 
required for such a simulation would be prohibitative, and so as a compromize we simply incorporate a gaussian 
energy input into the thermal history calculation, which models what we believe would be the effect of a network of 
active sources. The effective influence of a network is given by a smooth temperature change of the baryons over a 
certain period of time. A simple way to introduce a smooth “jump” in the temperature is by an errorfunction, so 
that the heating rate is a gaussian. This model based approach allows us to investigate whether there are potentially 
interesting effects, before resorting to the more time consuming simulation based approach. This source has three 
parameters, the redshift of the maximum energy input, the amplitude at the maximum and width of the gaussian 
which models the increase to and decrease from the maximum. Surprizingly, we find that the limits on the spectral 
distortions in the black-body spectrum of the CMB provided by the Far InfRA-red Spectrometer (FIRAS) instrument 
on the COsmic Background Explorer (COBE) satellite, only constrain these parameters very weakly. 

The thermal history calculation yields the so called photon visibility function Q, parameterized by time, which acts 
as a source for the linear Einstein-Boltzmann solver CMB FAST 134 . This function encodes statistical information 


on when the photons which we observe today were last scattered. For the standard thermal history including just 
recombination, it can be modelled as a gaussian centred around z ~ 1100 with width A z « 50. When we include 
the energy input there are four interesting situations which can occur. If the energy input is around or just after 
recombination then it is possible to modify the time and length of the last scattering epoch. If the energy input occurs 
once the recombination epoch is ostensibly over, and is sufficiently short for some of the photons to remain unscattered, 
then it is possible to have effectively two surfaces of last scattering, one around the time of recombination and the 
other around the time of reionization due to the energy input. If the energy input is sufficiently late for recombination 
to be complete, and long enough for almost all the photons observed today to be re-scattered at reionization, then 
it is possible for there to be a single surface of last scattering at a much lower redshift, effectively shifting the time 
of last scattering. Finally, if the period of heating is very late and long, then the universe remains at least partially 
ionized for most of the time after recombination and becomes totally ionized after z ~ 10. We shall describe these 
four possibilities as delayed, double and shifted last scattering, and total ionization in the rest of this paper. Of course 
none of these possibilities is totally fundamental and just about anything is possible for a sufficiently complicated 
source, but they do have some illustrative value. 

In the next section, we discuss the calculation of the thermal history. First, we include a detailed review of the 
standard thermal history used in most linear Einstein-Boltzmann solvers. Then we introduce our topological defect 
motivated energy source and illustrate the different effects it can have on the thermal history of the CMB by reference 
to the photon visibility function. The effect of this source on the black-body spectrum of the CMB is discussed 
and we show that the current limits on spectral distortions would have to be substantially improved before we could 
exclude such thermal histories. In section III, we discuss the effects of these modified thermal histories on the angular 
power spectra of temperature anisotropies and polarization, using simple analytic arguments to provide qualitative 
understanding and a linear Einstein-Boltzmann solver to give quantitative results. Finally, we discuss the possible 
implications for more realistic active models such as topological defects and causal white noise models. It should be 
noted that we have used natural units (Ti = k = c = 1) throughout this paper. 


II. THERMAL HISTORY CALCULATIONS 


A. The standard thermal history 


Originally, the thermal history of the CMB was studied using the Saha equation. This gave sensible quantitative 
results, which were subsequently extended by Peebles and Zel’dovich |l5],|6j] to include various corrections due to 
the complexity of recombination to the ground state of hydrogen. These calculations have been further extended to 
include more aspects of the underlying Boltzmann equations for the photons and electrons, and calculations are now 
at the stage where further improvements should only lead to about 1% corrections to the angular power spectra of 
temperature anisotropies and polarization j|] , although even further improvements [flTf lead to somewhat higher than 
1% corrections. In this section, we review these calculations of the standard thermal history. 

We quantify the change in the number density of a particular species in terms of the relevant Boltzmann equations. 
Strictly, speaking there are seven equations for the protons, hydrogen, helium, singly ionized helium, doubly ionized 
helium, electrons and photons. However, it is sufficient to treat the helium, both singly and doubly ionized, just 
using the Saha equation approximation, since the recombination rate for helium is much faster than the expansion 
rate during the relevant epoch 1381. Also we are really only interested in the evolution of the photons and electrons 
since the interaction of the photons with the relatively heavy baryons is negligible. The number of photons is much 
larger than the electrons and therefore their evolution can be decoupled from the electrons, with this interaction 
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being treated in terms of the spectral distortions discussed later. Hence, we can model recombination in terms of the 
fraction of ionized free electrons x e = n^ Be /n H , where n H = n B (1 — Y He ) is the number density of hydrogen nuclei, 
n B is the number density of baryons and Y He ss 0.24 is the mass fraction of primordial helium created at Big Bang 
Nucleosynthesis. In fact, 


x e = Zh + 


i 
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where ,th is the fraction of ionized hydrogen and xa e is that due to helium. Using the Saha equation, one can show 
that the fraction of He + is about 10 -5 for redshifts below z = 2000 and for He ++ this is even lower. Hence, iHe 
is always less than 10 -5 and can be included simply into the calculation using only hydrogen. One slightly odd 
side-effect of including helium in this way is that x e can be greater than one, although this is only the inclusion of 
helium modifying the calculation in the appropriate way; xr is always one or less. The equation which governs the 
recombination of protons and electrons into hydrogen is 


dx h 
dt 


= -r 1 +t~ 1 + t~ 1 

''rec “ ‘'pi ~ t 'ci 
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with tpj 1 and t ~ { 1 being the rates for recombination of hydrogen, photoionization and collisional ionization. In 
the standard case, there are no external sources of photoionization and collisional ionization and, therefore, we just 
need to calculate the recombination rate. 

However, to do this we will also need to model the evolution of the temperature of these distributions. The photon 
temperature is just redshifted by the expansion of the universe T 1 = T 0 (l + z ) at redshift z, where T 0 = 2.728K is 
the current temperature of the CMB, while all the other non-relativistic species remain in thermal equilibrium with 
each other. One can derive the evolution of the electron temperature T e using the first law of thermodynamics [fj~5f, 


dTe 

dt 


2 a Te+ 3 1- 31W4 + (1 - Y He ) x e ^ A) 
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where a is the FRW scale factor, T is the heating rate per baryon and A the cooling rate per baryon. The first term 
on the right hand side is cooling due to expansion, the third term characterizes cooling or heating due to the change 
in the number of free particles. The second term is a summation of the heating and cooling rates of the physical 
processes involved 


r — r src , a — A rec + a cmb , (4) 

where A rec is cooling due to recombination, A CMB is the Compton cooling from the interaction with the CMB photons 
and T src is any energy input which we might postulate, assumed to be zero for the standard case. This can be further 
simplified by realizing that the recombination cooling due to the loss of kinetic energy from the changing number of 
free particles is 


Arec — ^He)Uec ■ (3) 

In the standard case, where we have no source of photoionization or collisional ionization, one can use (Jsij) to replace 
t~ e l with dx e /dt , which then exactly cancels the last term in (|j)[. Physically, we are just cancelling off two processes 
which are in equilibrium. Hence, in order to calculate the standard thermal history we must just compute and 
A C mb, and then use a numerical routine to solve the differential equations for x e and T e , although the stiff nature of 
these differential equations does require some caution. 

In order to compute these two quantities we shall now assume a flat, Ho = 1, universe with the Hubble parameter 

H(a ) = H 0 a ~ 2 (a + a e J 1/2 , (6) 

where a 0 = 1, a eq = (l + ^ eq ) _1 and z eq = 2.40 x 10 4 /i 2 (T o / 2.728K) -4 is the redshift of radiation-matter equality, with 
the Hubble constant parameterized in the usual way, Hq = 100/ikmMpc -1 sec -1 . The recombination rate that we 


1 Strictly speaking the description in the text should include the effects of helium, but exactly the same cancellation will take 
place when this is included correctly for the same physical reasons. 
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must calculate here is the ‘net’ rate, taking in account both the recombination and ionization rate to and from all states 
of the hydrogen atom. Recombination directly to the ground state produces a Lyman-a photon, which, with high 
probability, immediately ionizes a hydrogen atom, either the one which it has come from or one of its close neighbours. 
Hence, we do not have to consider recombination to the ground state with the exception of two possibilities. Firstly, 
some of the Lyman-a photons may be redshifted out of their resonance line, before they have the chance to be 
reabsorbed. Failing that the ground state can only be reached by the two photon decay: 2si / 2 —> 2 pi/ 2 + 7 —> Is + 7 . 
The net recombination rate for transitions to and from states above the ground state and these two possibilities for 
the ground state is given by |35| , ^9|| 

t~l = a e x H n H C - P e (1 - %) e _3A/4Te C. (7) 


This complicated expression requires some explanation. Firstly, and most simply, we define A ss 13.60eV to be the 
binding energy of the ground state. The rate of recombination to all excited levels is ff(i|l 


a e = 2 A 




( 8 ) 


where A = 2 5 3 3 / 2 a 3 TrAg = 2.105 x 10 22 cm 2 given in terms of the fine structure constant a = 1/137 and the Bohr 
radius Ag = 0.529A, and m e = 0.511 MeV is the electron mass. Also included is a quantum correction/ for radiative 
effects known as the Gaunt factor which is given by g ~ 0.943 for temperatures below 5000K. The function <f>'(t e ) 
comes from summing up the interaction cross-sections of all excited levels and is given by ©13 

(te) = \ (^1.735 - Inte + |) - ^/^(l/ie) , (9) 

where 



( 10 ) 


is the exponential integral function and t e = T e / A. The ionization rate /3 e is related to the recombination rate by a 
detailed balance argument and local thermal equilibrium between all the excited states, and is given by 


fte — 



3/2 

e -B 2 /T e ) 


( 11 ) 


with B 2 = A/4 being the binding energy of the lowest lying excited, n = 2, state. Finally, the correction due the 
redshift of Lyman-a photons and the two photon decay is |35|,|39| 


1 + KVni s 
1 + K(V + P e )n u ’ 


( 12 ) 


where K = A 3 a/87ra, X a = 1216A is the wavelength of the Lyman-a photons, V = 8.23s _1 is the net rate of the 
two-photon decay and n i s ~ (1 — iH)n H the number density of hydrogen atoms in the Is state. At very late times the 
density of the baryons and electrons is low and direct recombination to the ground state is possible, since the density 
of the produced Lyman-a photons is then low as well [j42|. This leads effectively to C —> 1 at low redshifts, where we 
have choosen z < 100. The inclusion of this effect did not change the reionization histories we studied. 

In order to calculate the Compton cooling rate, one must use the Boltzmann equation for the photon distribution, 
which in full generality is jf^ [id]] 


dn . . dn 

TA A) = -di 
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~dt 
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br 


dn 

dt 


dc 


(13) 


where the subscripts refer to Compton scattering, production of bremsstrahlung (free-free) photons and double Comp¬ 
ton scattering, respectively and v is the frequency of the photons. We only need concern ourselves with the term due 


2 A more precise fitting formula for the recombination rate is given in ||. However, the effect on the anisotropy power spectrum, 
by including these corrections is less than 2%. 


4 











FIG. 1. An example of the energy input. On the left the energy input per free particle per comoving volume, q, as a function 
of redshift and on the right the corresponding rate of energy input, —dq/dz. 


to Compton scattering, since the other two processes lead to negligible cooling. When integrated over frequency, this 
gives the Compton cooling rate for Thomson scattering of hot electrons off photons in the plasma [^3|-^5|Jl2[ 
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where the Thomson scattering cross section is tr T = 6.65 x 10 29 m 2 and n c mb is the undistorted cosmic microwave 
background spectrum, n^ B = — 1. 


B. The gaussian energy input 


In the previous section we reviewed the standard thermal history of the CMB. We ignored the effects of photoioniza¬ 
tion, which is usually assumed to be the source of ionization in the intergalactic medium due to early object formation, 
and set to zero the heating due to sources, T src . As already stated the purpose of this paper is to investigate the 
possible effects of the active sources themselves and in this section we introduce a phenomenological expression for 
T src , which is intended to model the effects of a network of topological defects, specifically cosmic strings. In order to 
do this we have to make various assumptions which basically allow us to say that the whole universe becomes ionized 
in a essentially homogenous way. This is unlikely to be completely true in a realistic model, since sources are random, 
but it is required for us to make calculations possible. The effects of inhomogeneous reionization contribute only to 
second order in the CMB anisotropies and will be the subject of future work. 

We assume, therefore, that the active sources are distributed homogeneously in the universe and more importantly 
that the density of these sources is large enough that they will interact significantly with all the baryons over a short 
timescale. At early times the thermal velocity of the particles in the plasma is large and the velocity perturbations 
will be relatively small, so the sources will have negligible effect on the baryons. But as the perturbations grow and 
the thermal velocities are redshifted by expansion, the sources will become more significant. In the case of cosmic 
strings, it has been suggested 1331 that the formation of wakes starts around z ss 800, which can shock heat the 
plasma upto temperatures of a few million Kelvin during subsequent epochs, dependent on the small scale properties 
of the string network. The increasing kinetic energy of the plasma leads to a partial reionization of hydrogen, but 
as the density of active sources decreases, the energy release in the plasma will also decrease. Therefore, the heating 
takes place only over a finite period of time and afterwards the temperature of the baryons will remain constant, if, 
for the moment, we neglect the effects of cooling. One could model this effect in many ways and we have chosen to 
do this using a smooth function which interpolates between temperature of the plasma before the energy input and 
the temperature after heating by the sources, once again assuming no cooling. More specifically, all the baryons are 
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heated up to a temperature of Th ea t during a time interval of length A 2 = p, centred around z = 2 , using the energy 
per free particle per comoving volume 
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where erf(x) is the error function defined, in the standard way, by 


erf(x) = 



Therefore, the rate of energy input has a gaussian shape 


(15) 


(16) 


r 


src 



dq dz 
dz dt 


3 

2 




/-i . \5/2 -^O^heat 

(1 + 2 ) - 

P 


exp 




(17) 


where p is the width of the Gaussian, Theat/p is proportional to its height and 2 is the position of the peak. The 
energy input © and the rate (0) are plotted in Fig. 0 for 2 = 700, p = 200 and T heat = 5.0 x 10 6 K. We want to 
empasize again that our heating model is just an Ansatz, which is what we believe a reasonable one and allows one 
to predict the behaviour of the CMB caused by such a heating process. 

Of course, the plasma will never achieve these high temperatures since Compton cooling is very efficient for 2 > 10, 
and in fact we find that it is difficult for the actual electron temperature to get much above T e ss 5000 K. If the 
effects of the energy input are not significant at late times, then the temperature of the electrons will drop back to 
the CMB temperature, once the heating has stopped. However, it is possible for the electron temperature to remain 
around 5000K if the energy input is significant at low redshift (see, for example, the total ionization model discussed 
below). 

So we now have a source term to add to the standard thermal history which is motivated by a network of evolving 
topological defects. In the standard thermal history we made a number of assumptions, which need to be re-examined 
in the presence of this source. Firstly, there is still no source of photoionization, since ionization by free-free emissions 
is negligible for our shock motivated source. The addition of the source can also lead to collisional ionization, but we 
have explicitly checked that this is negligible. The main reason for this is that the densities are not high enough and 
the fact that A ss 13.60eV « 158000K is much higher than the temperature of the plasma achieved by the thermal 
input. Similarly we have checked that collisional ionization and excitation do not play a significant role in the cooling 
process since they are much weaker than Compton cooling. 

Hence, we conclude that to model the effects of a topological defect network on the thermal history of the CMB, one 
can just modify the standard calculation by the inclusion of the source term ( p~7| ) bearing in mind the uncertainties 
of the Ansatz. One might wonder how the ionization occurs, since we have explained above that the usual physical 
processes - photoionization and collisional ionization - are negligible. The basic mechanism is by modifying the rates 
for recombination and ionization, a e and /3 e respectively, creating a shift in the balance between atomic hydrogen 
and free protons and electrons. More specifically, the increase in T e creates a significant modification to a e (the 
probability that an electron is captured by a proton decreases with increasing T e ) and a slightly smaller effect in f3 e , 
reducing hence shifting the balance towards free protons and electrons. In other words, ionization dominates 
over recombination as long as the matter temperature is large. 

We have investigated the effects of using this source for a wide range of parameters and found that there are 
four cases which can illustrate interesting effects. We call them delayed , double and shifted last scattering, and total 
ionization and an example of each is discussed below. Before doing this we should discuss what we shall use to quantify 
their effect on the microwave background. The differential optical depth due to Thomson scattering is r = x e n H aTa 
and hence the optical depth of any particular conformal time is 


r{r.l) 



(18) 


where 770 is the conformal time of the present day. From these two quantities we can construct two photon visibility 
functions [Q, the last scattering visibility function 

g{V) = f(v)e- rM , (19) 
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FIG. 2. The effects of delayed last scattering on the thermal history and photon visibility functions. We have used z = 1000, 
p = 350 and Theat = 1.3 x 10 8 K. The plots are : (top,left) the electron temperature T e , (top,right) the fraction of ionized 
electrons x e , (bottom,left) the last scattering visibility function, (bottom,right) the cumulative visibility function, all plotted 
against redshift a. On each plot there are two curves, the modified (solid line) and standard (dotted line) thermal histories. 






FIG. 3. The effects of double last scattering using z = 500, p = 100 and Th ea t = 10' K. The plots are arranged in the same 
order as Fig. 2. 
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FIG. 5. The effects of total ionization using z = 350, p = 300 and Theat = 1.5 x 10 7 K. The plots are arranged in the same 
order as Fig. 2. 

















and the cumulative visibility function 


h{v) = e ~ T{ri) , 


( 20 ) 


which encode statistical information about when the CMB photons observed today were last scattered or influenced by 
objects, such as topological defects, along the line of sight. We have already commented that in the standard thermal 
history the last scattering visibility function is a gaussian of width A z ~ 50 centred around z ~ 1100. The cumulative 
visibility function is even simpler, being effectively zero for z > 1100 and increasing sharply to one for z < 1100 in 
the standard case. This reflects our inability to see directly any objects beyond the surface of last scattering since 
the universe is opaque. It also has a simple statistical interpretation; for any given conformal time rj, then h{rf) is the 
probability that a given photon was scattered before rj and 1 — h(rf) is the probability that it was scattered after 77 . 
We shall see that these functions occur naturally in the calculation of CMB anisotropies and hence it is important for 
us to understand the effects of our energy input on them. 

Delayed last scattering - If the energy input is significant during or close to the epoch of recombination then it is 
possible to delay standard recombination. An example of this is z = 1000, p = 350 and Th ea t = 1-3 x 10 8 K, whose 
effects on the thermal history and photon visibility functions are illustrated in Fig. [| We see that the temperature 
of the electrons deviates away from that of the photons around z « 1500 and only manages to get back down to be 
close to the photon temperature around z « 100. This allows a significant delay in recombination process with the 
ionization fraction remaining greater than 70% until around z « 700. The last scattering visibility function is still 
similar to a gaussian, but it peaks around z « 700, as opposed to z « 1100, with a larger width and consequentially 
smaller amplitude. A corresponding change in the cumulative visibility function is also observed, with the transition 
from zero to one happening around z « 700. 

Double last scattering - If the energy input is slightly later, then it is possible for recombination to take place in 
the standard way, with x e decreasing almost to zero before increasing once again as the electron temperature begins 
to rise due to the energy input. If the period of heating is short enough for a substantial fraction of the photons 
which were originally last scattered at the standard time to remain unscattered during the reionized epoch, then it 
is possible to have effectively two surfaces of last scattering. This can be achieved, for example, if z = 500, p = 100 
and Theat = 10 7 K, and it effects are compared to the standard case in Fig. [| As for delayed last scattering the 
electron temperature deviates significantly from that photon temperature, but now it takes place over a relatively 
short timescale, between z ~ 700 and z « 300. The first epoch of recombination is effectively over by around the 
time at which T e begins to increase, but this increase leads to ionization which reaches a maximum of around 30% by 
z « 450, before decreasing again. The last scattering visibility function can be approximated by two gaussians one 
centred around z ~ 1100 and the other around z ~ 500. The cumulative visibility function has two steps, from zero 
to 0.3, and then from 0.3 to one. This tells us that 30% of the photons observed today were last scattered around the 
time of standard recombination, while the other 70% were last scattered during the epoch of reionization. 

Shifted last scattering - If the energy input takes place after the standard recombination has taken place, as for double 
last scattering, but the duration for which it takes place is much longer, then most of the CMB photons will be 
re-scattered during the epoch of reionization and the photon visibility function will be almost zero around the time 
of standard recombination. Although, this is physically distinct from delayed last scattering, the photon visibility 
functions will be very similar apart from a shift to lower redshift, and hence we call this shifted last scattering. Using 
z = 500, p = 200 and T heat = 4 x 10 7 K, one can shift the surface of last scattering to lower redshift, as illustrated in 
Fig. [|. In this case, the electron temperature is above that of the photons from around z ~ 800 until very close to 
the present day (z ~ 10), with the fraction of free electrons being larger than 10% between z « 700 and z « 100. The 
last scattering visibility function is now approximately a single gaussian centred around z ~ 500, with a much larger 
width than in the standard case. The cumulative visibility function can be approximated by a single step at z« 500, 
re-enforcing the fact that almost all the photons were last scattered during the epoch of reionization. 

Total ionization - If the energy input is very late and over a substantial period then it is very difficult for the universe 
to become neutral to any great degree between the epoch of standard recombination and the present day. In this 
case, the photon visibility function will be relatively small since a large fraction of the CMB photons observed today 
have not stopped scattering through the whole history of the universe and effectively there is no well defined concept 
of last scattering. One interesting side effect of this is that the universe becomes totally ionized at around z « 10 
and the electron temperature does not come back down to the photon temperature, remaining at around 5000K. An 
example of this type of thermal history is given in Fig. I using z = 350, p = 300 and Ti iea t = 1.5 x 10 7 K. The electron 
temperature is much higher than the photon temperature for z < 700, where the energy input becomes significant, 
and remains around 3000K until around z ~ 20. At this point the effects of Compton cooling reduce substantially and 
the electron temperature increases to be greater than 5000K for z < 10. There is a corresponding two step increase 
in the fraction of free electrons. The last scattering visibility function is much smaller than in all the other cases, but 


9 


still has a peak around the time of standard recombination and a broader peak across a wide range of redshifts. The 
cumulative visibility function increases from zero to about 0.1 around the epoch of standard recombination and then 
increases very slowly to one between z « 500 and z ~ 1. Hence, only 10% of the photons observed today were last 
scattered during the epoch of standard recombination and the other 90% were scattered at some point after z ~ 500. 
The amount of energy per baryon injected into the universe is given by E/N\> = 3/2Th ea tj which gives in the case of 
delayed last scattering 16.8 keV per baryon. After the baryons reach equilibrium with the CMB this energy input is 
recognizeable as spectral distortion of the Planckian spectrum which is discussed in the next subsection. The resulting 
distortions should be small since there are about 10 9 photons per baryon. 

We should emphasize that none of these possibilities is fundamental and the effects described above will often be 
superposed in a non-trivial way. Also only the total ionization model achieves substantial ionization at late times and 
hence some other mechanism, probably photoionization, would be required to pass the Gunn-Peterson test. 

Although the mechanisms of heating and reionization discussed in this paper are entirely different, the results can 
look very similar to reionization due to the presence of light supersymmetric particles (e.g. a light photino or higgsino) 
|^8|,^9|. These light inos decay into UV-photons which subsequently ionize the matter. Dependent on the lifetime 
of light inos the (re-)ionization history is comparable to the one discussed in the present paper. The remaining 
energy of the ionizing photons can heat the matter to much larger temperatures as in the cases discussed here [ ff9| . 
Other mechanism of reionization which have been studied in the past are reionization by decaying massive neutrinos 
[fio| [52|| or the evaporation of primordial black holes |^,^|. It should be noted that these mechanisms can also lead to 
reionization histories similar to the ones presented in this paper. 


C. Spectral distortions 


We have already noted that the number of photons is so much larger than the number of baryons and hence the 
evolution of the photon spectrum can be decoupled from the evolution of the number density of free electrons during 
standard recombination and subsequent epochs of reionization. But the spectrum does evolve and, assuming the 
deviations from a pure black-body are small, as they will be for cases under consideration here, the evolution of the 
spectrum can be studied in terms of the spectral distortions, known as the Compton y-parameter or y-distortion 
[p|,[l6 54 5f|, the free-free distortion ]56| , |57[ and the /z-distortion | pi~i| . [T6| . [58| . The COBE satellite had on board an 
instrument, known as FIRAS, which measured the spectrum of the CMB to high degree of accuracy, placing apparently 
stringent upper bounds on these distortions. Subsequently, it was shown these upper bounds placed constraints on 
energy input at epochs before the standard time of recombination, but as we will show they place only very weak 
limits on the type of energy input that we have discussed in the previous section. 

The source of energy input we have introduced will only lead to free-free and y-distortions, with the y-distortion zero, 
if there is no significant energy input before the time of standard recombination, since the rate of Compton scattering 
is much slower than the expansion at these redshif ts [plfl . Similarly, we can ignore double Compton scattering relative 
to single Compton scattering after recombination |6l| l . Hence, we need only consider the y-distortion and the free-free 
distortion. Compton scattering conserves the total number of photons and hence the only way to increase the energy 
of the radiation due to the interaction with the hot plasma is by redistributing the photons from lower frequencies to 
higher frequencies. This is measured using the y-distortion. The thermal bremsstrahlung process on the other hand 
produces photons in the low frequency region, creating the free-free distortion. If the temperature is measured in the 
Rayleigh-Jeans part of the spectrum, it deviates due to the increased/decreased number of photons in this part of the 
spectrum by |i^ , |59|| 


Trj{x) - T 0 _ 

To 

Trj(x)-T 0 _ Y a 
To ~ x 2 


( 21 ) 


where Trj(x) is the Rayleigh-Jeans temperature and Tq is the measured temperature of the Planckian distribution. 

The COBE FIRAS limits |k| on these two parameters are |y| < 1.5x 10~ 5 (95%CL) and |Yj t | < 1.9x 10 -5 (95%CL). 
In fig. (j(|,left) we have plotted the brightness I (/ <x x 3 n(x), with x oc v/T) of a spectrum with y = 1.5 x 10 -2 . The 
redistribution of low frequency photons to higher frequencies, due to the interaction with the hot plasma, is evident 
from this graph. In fig. (|y, middle) the brightness of a spectrum with a distortion due to the release of bremsstrahlung 
photons from the plasma is shown, with Y tt = 0.1. Note that the distorted graph in fig. (||, middle) is only valid in 
the Rayleigh-Jeans limit. However, one can see the increasing brightness due to extra photons in the low frequency 
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region. Fig. (fright) shows the relative distortion of the brightness, I re \ = [I(x) — / P1 (a;)]// P1 (a;) where I pl (x) is the 
brightness of the Planckian spectrum. The values of the distortions in this plot are taken from the upper limits from 
the COBE FIRAS experiment. 





FIG. 6. On the left a y-distorted (solid line) and a Planckian (dotted line) spectrum are plotted, where it can be seen that 
the photons are redistributed from lower frequencies, to higher frequencies. In middle a free-free distorted spectrum (solid line) 
is plotted where the increased number of low-frequency photons is established. On the right the distortions with the value of 
the upper limits from the COBE FIRAS experiment are plotted (solid line y-distortion, dotted line free-free distortion) relative 
to the related Planckian spectrum ((n — npi)/npi)). 


The creation of spectral distortions is dictated by the photon Boltzmann equation ([[ 3 ]) sourced by the Kompaneets 
term for Compton scattering 


dn 

~dt 

and free-free bremsstrahlung via 


f T~ Id 

— n free cr„_ X _ 

6 m e x 2 dx 


T e dn 

Try 8X 


+ n( 1 + n) 


dn 

~dt 


[1 — n (e Xp - 1)] 


( 22 ) 


(23) 


br 


where x = 27tz//T 7 and x p = 2ni//T e are dimensionless frequencies relative to T 7 the unperturbed temperature of the 
CMB and the electron temperature T e . The coefficient k is given by, 


327 T 3 e 6 n ! ™ e n B g bt (x P ) 
3 m e T3 y/6-Km e T e 


(24) 


for a plasma containing both hydrogen and helium |d[] and g hr (x p ) is a Gaunt factor, accounting for quantum 
corrections to the free-free radiation process. One can either use a frequency averaged Gaunt factor g hI (x p ) = ~g hr « 1.2 
which gives an error less than 20% |i2|], or a more precise definition g br (x p ) = V3 hr (2.25 /x p ) /it for x p < 0.37 and 
9 hr{x p ) = 1 for x p > 0.37 fplf . 

To account for the distortions in the black-body spectrum, we solve the Boltzmann equation in terms of the 
unperturbed spectrum n q(x) and the spectral distortions y or Y {{ , that is, n{x) = no(x) + yf(x), where [H 


f(x) = 


-4 


(e x — l) 2 \tanh(a;/2) 

If we do this in an ft = 1 universe then the y distortion is given by 

_ v 

y = J 

to 


Tp — Tr 


m e 


Zh 

dt ft! 1.16 x 10~ n f l b h (1 - Y He ) J + z ) 1/2 dz , 


(25) 


(26) 


and the free-free distortion is given by 

th 

to 


‘'h z h _1/2 

Y a = J k Tb T Tl dt ft 2.33 x 10~ 6 g br flg/t 3 (1 - Yhs) J fl-^(l + z)^ x e (z){l + z) 1,2 dz 


( 27 ) 
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where ih and z h are the time and redshift when the energy input becomes significant. 

One can estimate an upper bound on the distortions due to our energy input, assuming that we achieve total 
ionization, that is, x e = 1 and a constant temperature, T e over a finite range, say z = Zh to z = 0, with z h >> 1. We 
find that the y-distortion is given by 0] 


y « 2.3 x 1CT 5 



and the free-free distortion is [0] 


Y« 


6.3 x 10” 12 



n b \ 2 / h \ 3 ( T e y 1 ' 2 / Zh y /2 
0.05 ) \0.5 ) ^ 0.76 yV 5000K / v 1000/ 


(28) 


(29) 


Since our heating model does not result in larger values than T e = 5000 K, it is easy to see that it would be difficult 
to achieve much more than y « 10~ 5 and Y {{ ss 10 - ' using standard cosmological parameters and therefore, the limits 
on spectral distortions from FIRAS do not constrain our model to a ny great degree. In fact, the actual spectral 
distortions are likely to be smaller than the upper bounds ( p8| ) and ( [29| ) since the ionization fraction and electron 
temperature are at their maxima only for a small portion of the range. We have calculated these distortions for the 
four models which we have been considering. They are y ss 3.2 x 10“ 6 and Y„ « 4.3 x 10“ 8 for delayed last scattering, 
y « 3.0 x 1(U 7 and Y t{ ss 5.2 x 10~ 9 for double last scattering, y « 1.7 x 10" 6 and Y ft « 2.6 x 10~ 8 for shifted last 
scattering, and hnally y « 9.8 x 10 -7 and Y ft « 1.6 x 10 -8 for total ionization. For comparable heating temperatures 
these are also the approximate distortions given in 00, however fi9| does not consider a particular model for the 
heat input. The heating temperatures in the model used in [ pT[ are higher than the ones discussed in this paper, 
resulting in much stronger limitations when reionization can occur. Similar results are obtained in 0f| under 
the assumption of stationary heating. Future measurements of the cosmic microwave background spectrum, like the 
experiments on the satellite missions MAP (Microwave Anisotropy Probe) 0], Planck Surveyor 0,0] and DIMES 
(Diffuse Microwave Emission Survey) 167 0, are at the level of 0.1% accuracy and can provide more stringent bounds 
on the distortions and could rule out some of the models discussed here. However this depends how accurately one 
can remove the galactic dust contamination from the data ||69| . 


III. ANGULAR POWER SPECTRA OF TEMPERATURE ANISOTROPIES AND POLARIZATION 

A. Simple analytic arguments 

To predict CMB anisotropies one has to solve the evolution equation for each species of particles either numerically 
or analytically 0 0]. One can understand the qualitative structure of temperature anisotropies and 
polarization using a formalism first applied to passive adiabatic models 00 and later adapted to incorporate 
active source models 0 0. For the purpose of this section we shall ignore the vector and tensor contributions to 
the anisotropies, and concentrate on the scalar component since it is this which is most affected by the modifications 
to the thermal history discussed earlier. In this case, the angular power spectra for the temperature anisotropies, Cp, 
and the polarization]^], Cp are given by 00 

q f‘ OO c\ nOO 

Cf = -/ k 2 dk(Ai(k,r] 0 )Ai{k,r ] 0 )), Cf = - / fc 2 dfc(^(fc, r,o)Ef (k, %)>, (30) 

n Jo n Jo 

where A(k,rj 0 , ^i) and E(k, r] 0 , /.i) are the temperature and polarization distribution functions expanded in terms of 
spherical harmonics of the angular variable y = cos 6 1 with coefficients Ai(k,rjo) and Ei(k,rjo). The basic procedure 
involves writing these multipoles in terms of an integral over the lower multipoles and the gravitational potentials 
along the line of sight. 

The line of sight integral for the temperature anisotropies is then [0 
rv o 

A/(fc, 770 ) = drj (S%(k, rj)ji [k (770 - v)} + S^k, r])ji° [k (y 0 - v)} + S%(k, [k (rj 0 - 77 )]) , (31) 

Jo 


3 Note that since we have ignored the vector and tensor components of the source, there is no magnetic component of the 
polarization 0. 
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where jf°(x) = j[(x), jf°(x) = (3 j"(x) + ji(x))/ 2, 

Sr(k, g) = g(g) (A 0 + ’F) + h(g) - $ ) , S^(k, g) = g(g)V B , <Sf (fc, g) = g(r))P , 

T, $ are the gravitational potentials characterizing the effects of the sources, V B is the baryon velocit 
(A 2 — Ei)/2. The equivalent expression for the electric component of the polarization distribution is |8(| 

rm 

E t (k,g 0 ) = -V6 / dgg(g)P(k,g)a [k{r] 0 - v)] , 

Jo 


(32) 


y^and P = 


(33) 


where 


3(l + 2)\j l (x) 

,,w = V 8 TTz? — 


(34) 


For the standard thermal history the last scattering visibility function can be approximated for pedagogical purposes 
by a delta function at 77 * and the cumulative visibility function by a step function at the same point, that is, 
< 7 ( 77 ) = 5{g — 77 *) and h(rj) = Q(g — 77 *). Hence, © and (|H|) can be approximated by 00 


A i(k, 770 ) = (A 0 + 'F) (k, g*)ji [k (vo - g*)} + V B (k, g*)ji [ k (go - g*)} + p (k, g *)Ji i k (go - g*)( 


+ f dr]U?-$\ ji [k (770 - 77 )] , 

J 77 * 


(35) 


and 


Ei(k,rj 0 ) = —V6P(g*) e i [k (go ~ g*)} ■ ( 36 ) 

Therefore, one can investigate the qualitative nature of the anisotropies and polarization by estimating Ao + H/, V B 
and P around the time of last scattering^] and T — <]> along the line of sight. 

We shall primarily be interested in the contribution from acoustic oscillations, since it is they which are most 
sensitive to the ionization history. Using the tight coupling approximation, that is, an expansion in powers of 1/r, 
one can deduce that {77J 


A 0 ( 77 *) = Aq(? 7 *) + $( 77 *) = e fc 2 /ks(»!.> 0 )Ao(0) cos ( AA. ) AAa 0 (0) sin 


kg* \ V3 - 

Vs) + ~‘ 

dr]'e- k2/k ‘ iri *' 71 ' ) sin ^-^( 77 * - 77 ')) F( 77 '). 


kg* 

Vs 


(37) 


where we have ignored the effects of baryons, F(g) = /c 2 ($ — H/)/3 is the structure function of the source, 


K 2 (t72,t?i) 


4 r 712 dg 

27 Jm VV)’ 


(38) 


is the damping length due to photon diffusion @00 and Ao(g) = Ao(g) + ^(g). This serves as an approximation 
for the intrinsic anisotropy created around the time of recombination. Furthermore, the other important quantities 

V B and P are related to A 0 by P ~ Vg ~ Ai ~ A 0 around this time. 

With a few subtleties, the three parts to the tight coupling solution ([u]) correspond to the passive adiabatic 

(A 0 (0) = 0,F(g ) = 0), passive isocurvature (A 0 (0) = 0,F(g) = 0) and active sources (A o (0) = 0, A o (0) = 0). When 
all the relevant effects are taken into account the contribution from Vb is suppressed relative to that from Ao + 4/, and 
that from P is even further suppressed. Therefore, there are in general two oscillatory components to the temperature 
anisotropy which are out of phase with each other, one which gives rise to peaks, Ao + T, and the other which fills 
in between them creating troughs, V B . These arguments lead to the much discussed acoustic peaks in the standard 
adiabatic scenario, which has the first peak at around l = 200 corresponding to the size of the sound horizon at time 


4 In fact, we will calculate the Ao + $ around the time of last scattering and assume that 4' — <f> is small. 
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of recombination. Isocurvature and active source models also have similar peak structure^], albeit with the main 
peak at slightly larger l |39|. The effective source of polarization Ai is out of phase with that for the peaks in the 
anisotropy spectrum and there is no contribution to the polarization from their source, Ao + T. This leads to a set 
of tight peaks which are out of phase with those for anisotropy. The amplitude of the polarization is generally much 
lower than the anisotropy since no net polarization is created during the tight coupling epoch |l],o,p4 


84 


We wish 

to modify this simple qualitative treatment of the structure of the anisotropies to the case where we have a more 
complicated ionization history. In order to do this, we replace the last scattering visibility function by 20 8 ^| 


g(r 1 ) = e TM 6(r] - 77 *) + (l - e tM ^J 5(r) - 7 ? r ), 


(39) 


where r/r is the time of reionization. Using this one can deduce that the cumulative visibility function will be given 
by h{rj) = 0 for 77 < 77 *, h{rf) = e~ T ^ for 77 * < 77 < r) r and h(rj) = 1 for 77 > r/ r . This is similar to the case of 
double last scattering discussed in the previous section which is the most general case. However, we shall discuss how 
this approach can be modified to understand the effects of the shifted and delayed last scattering scenarios. Since 
the visibility functions are not well represented in this way for the case of the total ionization scenario, we shall only 
comment briefly on its effects in this section. 

Using ([39]) we can deduce that [^0| 


A;(fc,?7o) = e 


— p ~ T (Vr) 


(Ao + '£) {k, 77 * )ji [k( 770 - 77 *)] + V B {k, 77 * )j t [k{r] 0 - 77 *)] + P(k, r] it )j l [k{r] 0 - 77 *)] 


+ (l - e~ TM ^j 

+ e~ rM J *7 - 6 ) ji [fc(? 70 -rf)]+ J d V “ $) 3l [Hvo ~ v)\ ■ 


(A 0 + '£) (fc, rj T )ji [k{rj 0 - 77 ,,)] + V B (k, r/ T )jl 0 [fc(?7o - r) T )] + P(k, 77 r )jf° [* 1(770 - 77 ,-)] 

Vo 

Vr 


(40) 


The first term in square brackets is due to photons which were last scattered at the time of recombination. Only 
a fraction e~ T ^ Vr ' > of the total number of photons observed today will have not scattered since that epoch and the 
remainder, a fraction 1 — e~ T ^ Vl \ was scattered at the time of reionization. The effect of these photons on the 
temperature anisotropy is given by the second term. The final two terms are a similarly modified version of the 
integrated Sachs-Wolfe (ISW) effect. 

We have already estimated the effects of Ao + 4/ and V B at the time of recombination, and these are just damped 
by the factor e _T ^ r - ) in these scenarios. However, we must now also estimate their effects at the time of reionization. 
We assume that the photons free stream from the time of recombination until just before reionization at time r] r — e, 
where e is usually small when compared to rj r . This yields 


A 0 (? 7 r - U k) = (A 0 + T) (k, 77 * )j 0 [k (r) T - e - 77 *)] + J dr 7 ('i' - j 0 [k (t] t - e - 77 )] . 


(41) 


If we now assume that the ISW component to this is negligible, as it will be in most applications, then this simple 
multiplication can be used as an initial condition for the epoch when the photons re-enter the phase where the photons 
couple again to the electrons. Using (p7|) . one can deduce that Q 


A 0 (77 r ) = A 0 (77 r ) + <E>(? 7 r) = 6 k * /k * 


A 0 ( 77 r - e) cos 


fee 

V3 


V3 y- , 

+ — A 0 (77 r - e) sin 


ke 

TV 1 


^ dri’e-b 2 /^^ sin 

^ «/r? r —e 


V3 


(?/r - 77') FW) ■ 


(42) 


Computing V B and P around the time of reionization is more tricky since there is pre-existing anisotropy. If this is 
large then it can modify the relation between V B: P and Ai. We do not believe that this will have a substantial effect 
on the temperature, but we shall return to this point when we discuss polarization below. 

We are now in a position to discuss the effects of this double last scattering scenario on the CMB anisotropies. 
The fact that there are two last scattering surfaces will lead to two sets of peaks with relative amplitudes e~ T ^ Vr ' > and 


5 We are only considering coherent active source models at this stage. 
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1 — e~ T(jlr \ whose scales will be set by the sound horizons at the time of recombination (A* ~ 77” 1 ) and the time 
of reionization. By substituting the free-streamed tight coupling solution into (| 42 |), one can deduce that this scale 
is k r ~ (77* + e) -1 , which implies that the size of the sound horizon is proportional to the time for which acoustic 
oscillations take place. Remember that e is small relative to ?7 r , but so long as 7y r 77*, then it is possible for e ~ 0 (v*), 
creating peaks which are based around very different scales. These scales can be projected into Z-space (Tt]] using the 
fact that the spherical Bessel function ji(x) is peaked around x = Z, which in the case of the anisotropy created at 
the time of last scattering is Z ss A* (770 — 77*) an d for that created at reionization is Z « Ay (770 — rj r ). 

This qualitative approach can be modified to also explain the nature of the structure of the anisotropy in the shifted 
last scattering scenario. In this case, the optical depth of the reionization is very large, r(?7 r ) 1 , which implies that 

one can ignore the first term in ([!(]), and hence the anisotropy is given by 

A i(A, Vo) = ( A 0 + ’F) (Vr )jl [Hvo - Vr)} + Vb{Vt )jl 0 [k(VO - Vr)} + P{Vr)jf° [HVO ~ Vr)} 

rm 

+ / dr] 

J T) t 

Furthermore, one can assume that the tight coupling regime effectively never ended at the time of recombination, 
which requires us to make the approximation rj T « e (and 77 * =0). This is clearly not totally true, since recombination 
of the protons and electrons did take place, but the only thing critical for estimating the anisotropies is the visibility. 
In the specific case of shifted last scattering which we are considering here, the visibility of the time of recombination 
is almost negligible. Therefore, we can estimate 


('F - ji [A (770 - 77 )]. 


(43) 


A 0 (77r) 


Ao(77 r ) + 4>(?? r ) = e -* 2 /*. a (%.0) 


A o (0) cos 






sin 


+ 0- 
(T ( , r _ n ,y 



(44) 


We see that now that there is just one set of peaks with their scale set by A r ~ t?” 1 . Since 77 r 77*, this corresponds to 

a shift of the entire peak structure to smaller Z. Also the diffusion damping length A“ 1 (?7 r , 0 ) has grown considerably 
and hence the effects of Silk damping j88| are prevalent on larger scales than in the standard scenario |^6) . 

We could also use this limit to understand the effects of delayed recombination, since there is just one surface of 
last scattering, but it seems more logical to treat this case as its name suggests as an increase in 77*, with r(?7 r ) = 0. 
This leads to a simple shift in the time of recombination, and hence the size of the sound horizon when the photons 
last scatter. The observational consequence is that the whole spectrum of anisotropies and also the Silk damping 
envelope are shifted to larger scales. 

The effects of these modified thermal histories on the position of peaks in the polarization is very similar to their 
effects on the temperature anisotropy. This is since the evolution of the source of polarization inside any surface of 
last scattering induced by reionization is very much the same as in the standard case. Therefore, if we make the same 
assumptions about the surface of last scattering, then the polarization is given by (2^] 

E\ ( k , 770) = —V6 je~ r ^P( 77 *)ez [£(770 - 77*)] + (l - e _T(l)r) ^ P{Vr)ei [k(vo ~ ? 7 r)]} • ( 45 ) 

As for the anisotropy, the function e; (x) oc ji (x) peaks around Z « x and the two contributions to ( | 45 | ) produce peaks 
on the two scales l « A*(770 — v*) and Z « k r (rjQ — rfr). Similar shifts in the peak position are possible in the cases of 
shifted and delayed last scattering. Since the polarization has only one source, as opposed to the temperature which 
has a number, this effect will be seen much more clearly. 

However, the amplitude of polarization on different scales is much more difficult to understand in any kind of 
generality. First, there are effects of Silk damping, which are likely to be very much the same as for the anisotropy. 
But more importantly there is the added difficulty of pre-existing polarization at the time of reionization. In the 
standard case, the quadrupole is negligible until very close to the time of last scattering and so its amplitude increases 
during that time, while oscillating out of phase with the monopole. It is the balance between this increase and the 
effects of Silk damping which makes the amplitude maximum at around the third or fourth peak. If the quadrupole is 
non zero before the time of reionization, as is likely to be the case particularly in the double last scattering scenario, 
then it is possible to shift this balance to larger scales. Hence, it is possible for the maximum amplitude to be at the 
first peak. 

It is hard to obtain any analytical predictions for the total ionization model. This is because the last scattering 
visibility function is smeared out over the whole history since recombination and can not be approximated by a simple 
function with two peaks (fig. ||). Therefore the line of sight integrals ([n]| 53 |) can not be carried out in a simple way. 
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The only prediction we can make is, since this scenario gives a small cumulative visibility up to very late times (see 
fig. (||)), last scattering visibility function), that the resulting anisotropy power spectrum will be suppressed up to 
relatively large scales. 




FIG. 7. The temperature (left) and polarization (right) anisotropy power spectra for an isocurvature CDM model. The 
dotted line refers to the standard thermal history, the long dashed line to delayed last scattering, the solid line to shifted last 
scattering, the dot-dashed line to double last scattering, and the short dashed line to the total ionization scenario. 


B. Isocurvature white noise 

Our main emphasis in this paper is on active source models. However, there is a class of passive models, known 
as isocurvature white-noise models, which have a number of the features of an active model. We have integrated the 
linearized Einstein-Boltzmann equations with an integrator known as CMBFAST |dJ for an isocurvature CDM model 
p0| p2[ with the initial power spectrum 

P t {k) = k\ (46) 

that is, an initial white-noise spectrum and the results for our sample of modified thermal histories discussed in the 
previous section are shown in fig. ( 0 )- The power spectrum is plotted up to l = 1500, because we do not expect 
interessting features on smaller scales in this first order approach. We should note first that, for the standard thermal 
history, the amplitude and spectrum of such a model is very much in conflict with the current observations [ p3| . 
However, as we shall discuss in section III.D, part of the motivation for this work is to investigate how the changes in 
the ionization history might create a more acceptable model. 

We see that the peak structure in the delayed and shifted last scattering scenarios is moved to much larger scales 
and that amplitude is suppressed by the damping envelope, very much in keeping with the analytic arguments of the 
previous section. For the case of double last scattering the analytic treatment suggests that there should be a second 
set of peaks at larger scales, but it appears that in this particular case that this second set of peaks is suppressed 
below that of the first set. The first peak which appears for the standard case at l « 350 also appears to be ‘washed’ 
out by the damping at the time of reionization. However, the three peaks which remain are on exactly the same scales 
as those in the standard case and represent the anisotropy created at the time of the initial epoch of recombination, 
suppressed by the damping factor e~ T ^ Vl \ Finally the total ionization scenario results in a heavy suppression of power 
up to very large scales (l ss 20) This is because the last scattering visibility function is smeared out over the whole 
time since recombination (fig.(j^)). We should note that reionization only appears to be effective on small scales and 
in each case the anisotropy on scales with l < 20 remains unchanged. 

As already discussed the peak structure of the polarization power spectrum is a more clean test of our analytic 
arguments, although our understanding of the amplitude is much less clear. In the standard case, the first peak is 
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around l « 250 and the damping becomes effective on scales with l > 1000. Once again the action of the delayed and 
shifted last scattering is easily understood from our analytic arguments with just a universal shift of the spectrum 
and damping scales. In the double last scattering scenario we now see that there are two clear sets of peaks which 
have different scales, those at the smaller scales being those created at the time of recombination damped by the 
appropriate amount. Since the visibility function at recombination is not completely zero, the damping phase leads to 
a quadrupole moment in the temperature anisotropies. This quadrupole contribution acts as a seed for the creation 
of polarization anisotropies ( |33| ) during reionization and hence the first peak in the polarization for this model is more 
prominent than in the others. Finally, total ionization results in a severely damped polarization spectrum, although 
even in this very extreme case one can just observe the remnants of the polarization created at the standard time of 
recombination. 


C. Simple coherent scaling source models 


In order to study active sources we will first discuss the simple coherent scaling source models introduced in 
refs. [^4 95 1 . Although these kind of sources are unlikely to be realized in the early universe, they have some some 
illustrative value. The scaling source is introduced as components of the Newtonian metric perturbations; specifically 
the curvature $ = <F 7 b + <F S and the gravitational potential \F = \F 7 b + \F S , where <F 7 b and >F 7 b are the contributions 
from the photon-baryon fluid, and <F S and \F S are those for the source. These sources terms can be related to the 
density of the source p s , its velocity v s and anisotropic stress n s by |^5| 


fc 2 $ s = 47rGa 2 + 3 ^v s /k^j , 
k 2 (*F S + <F S ) = —87rGa 2 7r s , 


(47) 


where the derivatives are with respect to conformal time p = f dt/a. The conservation equations for the seed source 
are 1951 


Ps + 3- (Ps+Ps) = ~kv s , 
a 

,a , 2 

v s + 4—u s = kp s - -fc7r s , 
a 3 


(48) 


where p s is the pressure of the source. Clearly, there are only two independent quantities here and we can in general 
choose any two arbitrarily, the other two being computed by evolving these equations. One simple scaling ansatz, 
which we shall call the pressure source, is given by |)5| 


4:TrGa 2 p s = p 1 ^ 2 


sin ( Akp ) 
{Akp) 


= 0 , 


(49) 


where 0 < A < 1 and another which we shall call the stress source has 


47rGa 2 7r s = p 1/2 

B'2 


6 


sin (Bikp) 
(Bikp) 


sin ( B 2 kp ) 
(■ B 2 kp ) 


(50) 


where 0 < (B i, B 2 ) < 1, with p s given by the pressure source. For the results presented here we have chosen to use 
A = 1.0, B 1 = 1.0 and B 2 = 0.5. 

Figs. (||) and (]|) show the results of incorporating the pressure and stress sources into CMBFAST. For the pressure 
source the anisotropy results follow the exact pattern predicted by our analytic arguments and already confirmed in 
the isocurvature white noise model. Now the main contribution to the intrinsic anisotropy is from the third term 
of the tight coupling solution (|37|), which moves the dominant peak in the temperature anisotropy to smaller scales. 
The delayed and shifted scenarios comprise a single set of peaks at larger scales than in the standard case, and in 
particular the main peak, which is at around l « 350 in the standard case, is moved to around l « 250 in the delayed 
scenario. In the case of double last scattering one can still recognize the peaks from the recombination epoch on small 
scales. However the expected second set of peaks from reionization is hardly recognizeable. The polarization in each 
of these cases also follows the pattern already established. 
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FIG. 8. Temperature (left) and polarization (right) anisotropy power spectra for a pressure source model with A = 1.0. The 
dotted line refers to the standard thermal history, the long dashed line to delayed last scattering, the solid line to shifted last 
scattering, the dot-dashed line to double last scattering, and the short dashed line to the total ionization scenario. 




FIG. 9. Temperature (left) and polarization (right) anisotropy power spectra for a stress source model with A = 1.0, B\ = 1.0 
and B 2 = 0.5. The key to the curves is the same as in fig.(^). 


One should note that in these truly active models the large-scale anisotropy is created by the ISW effect due 
to the existence of the sources along the line of sight. In the delayed, shifted and double last scattering scenarios 
reionization takes place very much before 2 : ss 100, whereas most of the ISW effect comes from sources present after 
this time. Therefore, reionization has very little impact on this contribution, that is, the penultimate term in ([h]) is 
effectively zero. In the total ionization scenario, where the universe is ionized for most of the time between standard 
last scattering and the present day, this is not necessarily the case and reionization can interfere with the anisotropy 
on large scales. This is illustrated in fig. |], although superficially it appears that the large-angle contribution has 
increased. This is not in fact the case since the normalization to COBE makes all the models equally around l = 10. 
What has happened is a redistribution of the large-scale power and a reduction of the contribution to the COBE 
normalization from scales smaller than l ss 10. This will not always be the case since it depends critically on the time 
when most of the anisotropy is created relative to the ionization history. 
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The results of using the stress source are very similar to those from the pressure source as illustrated in fig. |f]. 
However, there is one case which is quantitatively different. It appears that the first peak in the polarization spectrum 
for the shifted scenario is higher than that for the double last scattering scenario, where in the two cases already 
considered (the isocurvature white noise and pressure models) these peaks have almost the same height. We believe 
that this subtle effect is due to the quadrapole of the temperature anisotropy, which is the source for polarization, 
being non-zero when the second tight-coupling epoch begins, although we have no analytic reasoning for this. 


D. Causal white noise 

One of the original motivations for studying these non-standard thermal histories was to attempt to rectify some of 
the observational problems of a class of structure formation models known as causal white noise (CWN) models 
These models were spawned out of the realization that standard scaling models with defect motivated stress-energy 
components are unable to explain the observed matter fluctuations on 100h -1 Mpc scales in a Einstein-de-Sitter 
universe, unless in case where large scale biases are acceptable. It was suggested that if the source was switched off at 
some point before a critical redshift (z c ~ 100), then the power on 100h -1 Mpc scales was exactly that observed and 
the excess of power on smaller scales can be rectified by modifications to cosmological parameters. An observationally 
unacceptable side-effect of this is the spectrum of CMB anisotropies produced. 



FIG. 10. Temperature anisotropy power spectrum for a causal white noise model. The dotted line is for a standard thermal 
history, the solid line for a total ionization history with the parameters z = 400, p = 350 and The at = 3 x 10 7 K and the long 
dashed line is for the same model with inclusion of a cosmological constant Ha = 0.65, fib = 0.1 and f2 C = 0.25. The spectra 
include scalar, vector and tensor contributions and are normalized to COBE. Also included are the observations data points 



Fig. [h] illustrates these problems for a simple CWN model in the standard thermal history, plus the effects of our 
total ionization scenario, in both an Einstein-de-Sitter cosmology and also one with a non-zero cosmological constant. 
In the standard scenario we see that there appears to be an excess of power on small scales and the shape of the 
spectrum on large scales is in conflict with the COBE data. This is because these models effectively tilt the spectrum 
toward smaller scales. The action of reionization is to reduce the power on small scales to a more acceptable level, but 
the problems on large-scales remain. These can be partially relieved by the inclusion of a cosmological constant, but 
at the expense of reducing the power on smaller scales. It was concluded in ref. |j(],|97j that although such models can 
explain the formation of structure, they have considerable problems explaining the observed CMB power spectrum. 
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FIG. 11. The temperature anisotropy power spectra for a cosmic string model under inclusion of a cosmological constant 
Qa = 0.80. The key of the plot is the same as in fig. [t| We have also included the data points with errorbars |98[. 


E. Cosmic strings with a cosmological constant 

The other motivation is to investigate whether acoustic peaks in defect models may shift to larger angular scales 
in realistic thermal histories. We have already mentioned that these models appear to predict a peak (if any at all, 
see refs. |m| |26||io|] ) on much smaller scales than in the standard adiabatic scenario, and in fact only very convoluted 
models can rectify this [ pd| . One might ask is this a generic phenomena and clearly our earlier arguments suggest 
that the ionization history can be modified to allow this to happen in more generic defect models. 

We have applied the same modified thermal histories to the model of structure formation by cosmic strings with 
a cosmological constant presented in ref. |j29f . The introduction of a non-zero cosmological constant can improve 
the shape of the matter power spectrum and its amplitude on scales of 100h -1 Mpc to an acceptable level, and it 
was shown that an important consequence of this is a broad peak in the CMB power spectrum on around l = 500. 
The model presented here has fie = 0.05, fl c = 0.15 and = 0.80. In fig. [Tl| we have plotted the temperature 
anisotropy power spectrum for our thermal histories with this model. The spectrum behaves very much as one might 
have expected from our earlier analytic arguments and also like the simple coherent models, although it should be 
noted that in this incoherent model the concept of peaks in slightly different. For standard recombination (dotted 
line) the peak is at l ss 500 and has an amplitude of « 95 /iK, which has been shifted, for example, in the case of 
double last scattering (dot-dashed line) to a broader peak at l ~ 200 and a height of 60 /xK. Similar modifications are 
made in the shifted and delayed last scattering scenarios. This is clearly an improvement with respect to the current 
observational data, although maybe a not a necessary one. 


IV. DISCUSSION AND CONCLUSION 


We have introduced a heating source motivated by the active character of structure formation in the context of 
sources like a cosmic strings. The main assumption is that there exists a phase in the universe after recombination, 
where the density of active sources is large enough to heat the baryons homogenously, at least over a short period 
of time. After this the influence of the active sources on the heating was assumed to be irrelevant for our effective 
description. Sources of ‘late’ reionization, such as, photoionization through early object formation are not included in 
our analysis. We also have not considerd the second order contribution from the Vishniac effect 


100], which will 


lead to an extra contribution in the temperature anisotropy power spe ctrum at large l (l > 2000). An inhomogenous 
treatment would also have led to such second order contributions [ 1011. 

We calculated the spectral distortions for a wide range of heat source parameters and established that it is hard 
to violate the COBE FIRAS limits on these quantities. We extended the framework of the semi-analytical Hu and 
Sugiyama formalism ||77| to the case of reionization and analyzed four generic heat source types explaining their 
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influence on the CMB temperature anisotropy and polarization power spectra for isocurvature white noise and a 
simple scaling source models. The behaviour of the models was explained well by these arguments. We have found 
that if there is only a shift of the surface of last scattering, the anisotropy power spectrum becomes damped and 
the peaks are shifted to larger scales. If there appears a second surface of last scattering the acoustic peak structure 
changes and the suppression on small scales is not as large as for a just shifted last scattering surface. The most 
important feature of the polarization power spectrum is the appearance of a prominent contribution on intermediate 
scales due to reionization. We then applied the source models to CWN models and also to a realistic cosmic string 
model. We found that it was possible to reduce the amount of power on small scales in CWN models and to move the 
peak in string models to larger scales — the original motivation for this work. However, it appears to not be possible 
to have a substantial effect on the spectrum at large angles in CWN models. 
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